Marriage a-la-MOND: Baryonic dark matter in galaxy clusters and 

the cooling flow puzzle 

Mordehai Milgrom 

Center for Astrophysics, Weizmann Institute 



Abstract 

I start with a brief introduction to MOND phenomenology and its possible roots in cosmology-a notion that may 
turn out to be the most far reaching aspect of MOND. Next I discuss the implications of MOND for the dark matter 
(DM) doctrine: MOND's successes imply that baryons determine everything. For DM this would mean that the puny 
tail of leftover baryons in galaxies wags the hefty DM dog. This has to occur in many intricate ways, and despite 
the haphazard construction history of galaxies-a very tall order. I then concentrate on galaxy clusters in light of 
MOND, which still requires some yet undetected cluster dark matter, presumably in some baryonic form (CBDM). 
This CBDM might contribute to the heating of the x-ray emitting gas and thus alleviate the cooling-flow puzzle. 
MOND, qua theory of dynamics, does not directly enter the microphysics of the gas; however, it does force a new 
outlook on the role of DM in shaping the cluster gasdynamics: MOND tells us that the cluster DM is not cold dark 
matter, is not so abundant, and is not expected in galaxies; it is thus not subject to constraints on baryonic DM in 
galaxies. The mass in CBDM required in a whole cluster is, typically, similar to that in hot gas, but is rather more 
centrally concentrated, totally dominating the core. The CBDM contribution to the baryon budget in the universe is 
thus small. Its properties, deduced for isolated clusters, are consistent with the observations of the "bullet cluster" . 
Its kinetic-energy reservoir is much larger than that of the hot gas in the core, and would suffice to keep the gas hot 
for many cooling times. Heating can be effected in various ways depending on the exact nature of the CBDM, from 
very massive black holes to cool, compact gas clouds. 
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1. introduction 

Recent observations of the cores of galaxy clusters 
have undermined the long standing "cooling flow" 
paradigm. These observations do not detect the ex- 
pected telltale signs of cooling, and thus point to 
some, still moot, mechanism that heats the core gas 
and keeps it at an elevated temperature despite its 
short cooling time (for recent reviews see, e.g., Bauer 
ct al. 2005; Peterson and Fabian 2006; Sanderson 
ct al. 2006). Several heating mechanism have been 
proposed. The most popular, at present, seems to 
be heating by central AGN activity (See, e.g., Peter- 
son et al. 2003, Omma et al. 2004, and Nipoti and 



Binncy 2005 for a discussion of its pros and cons). 
Dark matter (DM), in its presently favored variety 
of weakly interacting particles, is not deemed a fac- 
tor in these considerations, as it does not affect the 
gasdynamics in the core (but see, e.g., Chuzhoy and 
Nusser 2006 for possible effects of more strongly in- 
teracting DM). 

MOND was introduced as an alternative to New- 
tonian dynamics that seeks to explain the dynamics 
in galactic systems without DM. It, arguably, works 
very well for galaxies, galaxy groups, and super clus- 
ters, and has made considerable theoretical headway 
(for reviews see, e.g., Sanders and McGaugh 2002, 
Scarpa 2006, Bekenstein 2006, Milgrom 2008). But 
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in galaxy clusters MOND does not completely ex- 
plain away the mass discrepancy. This is particu- 
larly so in the very cores of clusters where a large re- 
maining discrepancy is found. So even with MOND 
one has to invoke yet undetected matter peculiar to 
clusters: presumably in some baryonic form, here- 
after CBDM. Massive neutrinos with mass near the 
present experimental upper limit of 2 ev were also 
proposed (Sanders 2003, 2007). 

MOND determines the overall potential field in 
the cluster, but does not directly affect the micro- 
physics of the cooling flow and core gas. It does, 
however, shed new light on the possible role of DM 
in cluster gasdynamics: (i) Unlike the standard view 
according to which the DM constituents are weakly 
interacting particles, in the context of MOND the 
cluster DM is likely to be baryonic; so, it can possi- 
bly interact with the gas in ways relevant to the cool- 
ing puzzle, (ii) According to MOND, this CBDM is 
peculiar to clusters; so, it is not subject to the strin- 
gent constraints we have on ubiquitous baryonic DM 
candidates in galaxies (e.g., Carr and Sakcllariadou 
1999, Carr 2000, Kamaya and Silk 2002). (iii) Al- 
though in the inner parts of clusters, the required 
amount of CBDM is large compared with that ob- 
served in stars and gas, the total amount required 
for the cluster at large is comparable to the mass 
of the hot gas. This renders the baryonic option 
more palatable. It also greatly relaxes constraints 
on baryonic DM in clusters based on standard dy- 
namics, which requires 5-10 times more DM. Such 
small quantities of a new baryonic component easily 
satisfy nucleosynthesis constraints. 

I discuss the possibility that such a baryon com- 
ponent could keep the core gas from cooling precipi- 
tously, thus helping alleviate the "cooling flow" puz- 
zle. The total kinetic energy budget of the CBDM is 
known, and is sufficient to supply heating for many 
Hubble times. The idea of heating the cluster gas 
with baryonic DM is not new (e.g. Walker 1999), but 
as explained above, it takes quite a different shape 
in the context of MOND. 

I start by reviewing in section 2 the basic phe- 
nomenology of MOND, and its possible origin in 
cosmology; its significance and implications for the 
DM paradigm are discussed in section 3. In section 
4 I summarize what MOND tells us about DM in 
clusters. In section 5 I list candidates for the CBDM 
and discuss heating mechanisms. 



2. An overview of MOND phenomenology 

MOND can be described as a modification of dy- 
namics at low accelerations (Milgrom 1983a). More 
precisely, the MOND paradigm rests on the follow- 
ing premises: (i) There appears a new constant in 
dynamics, a , with the dimensions of acceleration, 
(ii) In the formal limit clq — > classical (pre-MOND) 
dynamics is restored, (iii) For purely gravitational 
systems, in the deep-MOND limit, a — > oo, the lim- 
iting equations can be brought to a form in which 
the constants do and G, and all masses in the prob- 
lem, rrii , appear only as rriiGao (Milgrom 2005) . This 
is a fiat that guarantees the required MOND phe- 
nomenology in purely gravitational systems. 

These requirements can be incorporated in vari- 
ous MOND theories. For example, in the nonrela- 
tivistic regime one can modify the Poisson equation 
for the gravitational field (Bckenstein and Milgrom 
1984), or one can modify the kinetic action of parti- 
cles leading to modified inertia (Milgrom 1994a). It 
follows from the basic assumptions that any MOND 
theory must be nonlinear even in the nonrelativis- 
tic regime (see e.g., Milgrom 2008 for details). From 
the above MOND tenets (and the assumption that 
ao is the only new dimensioned constant) also follow 
some important scaling laws (under which Gclq is 
invariant) for the deep MOND limit of purely gravi- 
tational systems in any MOND theory (see Milgrom 
2008 for more details): The theory is invariant under 
scaling of the radii r — > Ar, and times t — > Xt (masses 
are unsealed). (This is part of the conformal invari- 
ance, in space alone, of the deep MOND limit in a 
particular formulation of MOND-Milgrom 1997). 

As a consequences of this scaling invariance we 
deduce that if r(t) is a trajectory of a body in a 
configuration of masses to, at positions rj(t), then 
r(t) = Xr(t/X) is a trajectory for the configuration 
where are at Ar^i/A), and the velocities on that 
trajectory are v(i) = v(t/A). There is another scal- 
ing property: if we scale all the masses, all trajectory 
paths remain the same but the body traverses them 
with all velocities scaling as m 1 / 4 . 

Ultimately, one would like to incorporate these 
principles in a relativistic extension of General Rel- 
ativity. The state of the art of this effort is the 
TeVeS theory and its derivatives (Bekenstein 2004, 
and Bruneton and Esposito-Farese 2007). I shall not 
discuss theories here; for a recent review see Beken- 
stein (2006). The theories proposed so far involve 
a free function that essentially interpolates between 
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the above two limits. 

Schematically, ao plays a similar role to c in rel- 
ativity, or to h in quantum theory. On one hand, 
they all mark the boundaries between the classical 
and the modified regimes; so formally pushing these 
boundaries to the appropriate limits (c — > oo, h — > 

0. or ao — * 0) one restores the corresponding classi- 
cal theory. On the other hand, these constants also 
feature prominently in the physics of the modified 
regime and appear in various phenomenological re- 
lations. Examples of such relations for h are: the 
black body spectrum, the photoelectric effect, the 
Hydrogen atomic spectrum, the quantum Hall ef- 
fect, etc.. For c we could mention the (relativistic) 
Doppler effect, the mass-velocity relation, the life- 
time-velocity relation, and the radius of a black hole. 
These relations, in themselves, are independent in 
the sense that they do not follow from each other 
and are made related only through the underlying 
theory. The MOND paradigm similarly predicts a 
number of laws relating to galactic motions, some of 
which arc qualitative, but many of which are quan- 
titative and involve ao; they may be likened to Ke- 
pler's laws of planetary motionqjj. More details of 
these predictions and how they follow in MOND can 
be found in Milgrom (2008). Here is a partial list of 
such laws. 

1. The rotation curves for an isolated mass is asymp- 
totically flat: V(r) -»• (Milgrom 1983b). This 
follows for all MOND theories from the above 
mentioned scaling of length and time. 

2. The mass- velocity relation (aka the baryonic TF 
relation) between the total mass of a body and 
its Voo-.V^ = MGa (Milgrom 1983b, see anal- 
ysis in McGaugh 2005a, and Fig. [I] here). This 
also follows in all MOND theories from the above 
mentioned scaling of orbital velocities with mass 
in the deep MOND limit. 

3. In disc galaxies with high central accelerations the 
mass discrepancy appears always around the ra- 



1 The nominal Kepler laws were all discovered before the 
underlying theory of Newtonian dynamics was propounded. 
From this theory additional "Kepler's laws" can be deduced, 
such as the dependence of the Kepler constant on the mass 
of the central star. MOND was constructed by assuming 
only one law-the asymptotic flatness of rotation curves- 
based on rather anecdotal evidence existing at the time, 
and it was helped by a rough constraint from the Tully- 
Fisher relation. MOND then elevated these two into absolute 
predictions (replacing the original, velocity-luminosity Tully- 
Fisher relation by one between velocity and total mass, and 
predicting the exact power in this relation, still moot at the 
time). 
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Fig. 1. Galaxy mass vs. asymptotic velocity. Left: analog to 
the traditional Tully-Fisher plot involving stellar mass alone. 
Right: gas mass is included. The solid line has the log-log 
slope of 4, predicted by MOND and is not a fit (McGaugh 
2005b). 

dius where V 2 /R = gq. In galaxies whose central 
acceleration is below ao (low surface brightness 
galaxies) there is a discrepancy at all radii (Mil- 
grom 1983b, for tests see Sanders and McGaugh 
2002, McGaugh 2006). 

4. For a concentrated mass, M, well within its tran- 
sition radius, r t = (MG/ao) 1 ^ 2 , r t plays a special 
role (somewhat akin to that of the Schwarzschild 
radius in General Relativity). For example, a shell 
of phantom DM may appear in the vicinity of rt 
(Milgrom and Sanders 2007). 

5. Isothermal spheres have mean surface densities 
£ < a /G (Milgrom 1984) underlying the Fish 
law for quasi-isothermal stellar systems (see dis- 
cussion in Sanders and McGaugh 2002). 

6. A mass- velocity-dispersion relation <r 4 ~ MGcio 
underlying the Faber- Jackson relation for ellipti- 
cal galaxies (Milgrom 1984, 1994b). 

7. There is a difference in the stability properties of 
discs with mean surface density £ < ao/G and 
with a Q /G < E (Milgrom 1989a, Brada and Mil- 
grom 1999b, Tiret and Combes 2007). 

8. The excess acceleration that MOND produces 
over Newtonian dynamics, for a given mass, can- 
not much exceed ao (Brada and Milgrom 1999a, 
as confirmed for a sample of disc galaxies by 
Milgrom and Sanders 2005). 

9. An external acceleration field, g e , affects the in- 
ternal dynamics of a system imbedded in it. For 
example, if the system's intrinsic acceleration is 
smaller than g e , and both are smaller than ao, 
the internal dynamics are quasi-Newtonian with 
an effective gravitational constant ss Gao/g e 
(Milgrom 1983a, 1986a, Bekenstein and Milgrom 
1984, with applications in e.g., Brada and Mil- 
grom 2000a, b, Angus and McGaugh 2007, and 
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Wu et al. 2007). Some aspects of this prediction 
violate the strong equivalence principle, and thus 
conflict with the predictions of DM. 

10. The thin lens approximation breaks down in 
MOND, which also conflicts with the prediction 
of DM. 

11. Disc galaxies have a disc "DM" component in 
addition to a spheroidal one (Milgrom 2001). See 
confirming analysis in Milgrom(2001), Kalberla 
et al. (2007), and Sanchez-Salcedo, Saha, and 
Narayan (2007). This prediction conflicts with 
the expectation from CDM. 

12. Negative density of "dark matter" is predicted 
in some locations (Milgrom 1986b). Again this 
conflicts with predictions of DM. 

These predictions are all independent in the sense 
that one can construct galaxy models of baryons 
plus DM such that any subset of these relations are 
satisfied but not any of the others. Thus, in the con- 
text of DM they will each require its own separate 
explanation. In fact, as I indicated, some of these 
have implications that conflict with the predictions 
of CDM, and among them some that conflict with 
any form of DM. 

But, the flagship of MOND phenomenology is the 
detailed prediction of the full rotation curves of in- 
dividual disc galaxies from their observed baryonic 
mass alone. Its importance was evident from the 
outset (Milgrom 1983b); but, testing this prediction 
had had to await the advent of extended rotation 
curves afforded by HI observations: Rotation curve 
analysis in MOND started only some five years af- 
ter it appeared, with Kent (1987) and the rectifying 
sequel by Milgrom (1988). Some of the subsequent 
studies where by Begeman et al. (1991), Sanders 
(1996), Sanders and Verheijen (1998), de Blok and 
McGaugh (1998), Bottema et al. (2002), Gentile et 
al. (2004), Corbelli and Salucci (2007), Sanders and 
Noordermeer (2007), and Barnes et al. (2007). There 
are now of the order of 100 galaxies for which MOND 
has been tested in this way: one uses the observed 
baryonic mass in a disc galaxy to predict the ob- 
served rotation curve according to MOND (convert- 
ing stellar light to mass using the M / L value-mass 
to light ratio-a free parameter). It is, of course, the 
complete failure of this procedure with Newtonian 
dynamics that leads to DM in galaxies. MOND is 
performing extremely well in this regard. I show ex- 
amples in Fig. [5] for three galaxies of rather differ- 
ent properties: from the low mass, low acceleration, 
gas dominant NGC 1560, through the intermedi- 
ate NGC 3657, to the high mass, high acceleration, 



stellar-mass-dominated NGC 2903. More examples 
are shown in Fig. [3] Not only are the full rotation 
curves well fit with the above single parameter, but 
in gas dominated galaxies even this parameter is al- 
most immaterial, and MOND then practically pre- 
dicts the rotation curve. When the M/L parameter 
is important, its values, obtained with the MOND 1- 
parameter fits, are in good agreement with theoret- 
ical, population-synthesis estimates (e.g., Sanders 
and Verheijen 1998). 



2.1. a Q ? 

The constant ao appears in several of the above 
MOND laws of galactic motions, and its value 
was initially determined, with consistent results, 
by appealing to them (Milgrom 1983b). However, 
the best leverage on clq comes now from rota- 
tion curve analysis. For example, the value found 
in Milgrom (1988) was a sa 1.3 x 10~ 8 cm s -2 , 
and Begeman et al. (1991) found, with better 
data, a ps 1.2 x 10~ 8 cm s~ 2 (both with H = 
75 km s _1 Mpc~ 1 ). It was noted in Milgrom (1983a) 
that ao ~ and because of the still mys- 

terious cosmic coincidence (Qa ~ 1) we also have 
a « c(A/3) 1 / 2 /2^. 

In the context of quantum theory, the Planck 
length and the Planck mass, constructed from h, G 
and c, tell us where we can expect combined effects 
of strong gravity and quantum physics. Similarly, 
ao defines a length scale, £o = c 2 /a ps 10 29 cm, 
and a mass, M = c 4 /Ga ps 6 x 1O 23 M , that 
tell us where to expect MOND effects combined 
with strong gravity. In view of the above coinci- 
dences of ao with cosmological acceleration param- 
eters, £o and Mq are of the order of the Hubble 
radius £ « 2tt£ h , with £ H = cH^ 1 ps c(A/3)- 1 / 2 , 
and the cosmological mass parameter ("the mass 
within the horizon"), Mq ps 2ttMjj, with Mjj = 
c 3 G*- 1 (A/3)- 1 / 2 ps cPG^Hq 1 , respectively. This 
means that combined effects are expected only for 
the universe at large, and that there are no strongly 
MONDish black holes. 

These coincidences may point to a connection be- 
tween MOND and cosmology: either the same new 
physics explain MOND and the "dark energy" ef- 
fects, or MOND is an expression in local physics of 
whatever agent is responsible for the observed cos- 
mic acceleration (as embodied in A). If the param- 
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Fig. 2. The observed and MOND rotation curves (in solid lines) for NGC 3657 (left), NGC 1560 (enter), and NGC 2903 (right). 
The first from Sanders (2006), the last two from Sanders and McGaugh (2002). Points are data, dashed and dotted lines for 
the last two galaxies are the Newtonian curves calculated for the stars and gas alone; the reverse for the first (they add in 
quadrature to give the full Newtonian curve). 




Fig. 3. Additional MOND rotation curves from Sanders (1996) and de Blok & McGaugh (1998) (left) and from Sanders & 
Verheijen (1998) (right). (MOND curves in solid; stellar disc Newtonian curves in dotted; gas in dot-dash; and stellar bulge in 
long dashed.) 



ctcr a c = Mjj/£ 2 H = cH w c(A/3) 1 / 2 is somehow 
felt by local physics, it may not be surprising that 
dynamics is different for acceleration above and be- 
low this value. But why is it that cosmology enters 
local physics through its characteristic acceleration 
and not, for example, through its characteristic mass 
or length? I discussed possible reasons in Milgrom 
(1999). Note that a similar thing happens in the way 
earth's gravity enters the dynamics in experiments 



close to its surface: There, a "constant of nature"- 
the Galilei free-fall acceleration, g = M^G/Rm, ap- 
pears in the dynamics. If we know the escape speed, 
v es = (2M®G / '-Re) 1 / 2 , or the orbital speed on a 
near earth orbit v or b = (MqG / R^y^-which can 
serve as possible analogs of c^we would note the re- 
lation between this, g, and which is essentially 
the same as that between, c, an, and In. We would 
then probably deduce that the relation is a sign that 
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Galilean free fall is an effective law resulting from a 
deeper theory, which, of course, it is. 

Its relation with cosmological parameters raises 
the possibility that ao varies with cosmic time (Mil- 
grom 1989b, 2002). If it is always related to cHq in 
the same way, or if it is related to the density of "dark 
energy" , and this changes, then ao would follow suit. 
It is also possible, however, that A is a constant and 
so is ao- The possibility of a variable oo opens up in- 
teresting possibilities. Such variations may, in turn, 
induce secular evolution in galactic systems; for ex- 
ample, the mass-velocity relations dictate that the 
velocities in a system of a given mass should decrease 
with the cosmological decrease in ao (like a^ 4 in 
the deep MOND regime) , and this would be accom- 
panied by changes in radius, since adiabatic invari- 
ants, such as possibly rv, have to stay fixed (Mil- 
grom 1989b). Such variations in ao could also be in 
the basis of the cosmological coincidence Cl\ <~ 1 via 
anthropic considerations (Milgrom 1989b, Sanders 
2001). Also, the general connection of MOND with 
cosmology, and, in particular, the possibility of a 
varying, may provide a mechanism for inducing a 
local arrow of time by the cosmological one. For ex- 
ample, if in an expanding universe ao decreases, this 
fact is felt by small local system, and can induce a 
preferred direction of evolution in them. 

Analyzing the data of Genzel et al. (2006) on the 
rotation curve of a galaxy at z = 2.38, I find that 
they are consistent with MOND with the local value 
of ao; but, the large error margins, and the fact 
that at the last measured point the galaxy is only 
marginally MONDish, still allow appreciable varia- 
tions of a with red-shift. 

MOND, as it is formulated at present, does not, 
in my opinion, provide a clear-cut tool for treating 
cosmology (including the appearance of cosmolog- 
ical DM). It is possible that the understanding of 
cosmology within MOND will come together with 
the understanding of MOND within cosmology, in 
the sense alluded to above. 

3. Significance for the DM paradigm 

Can the successes of MOND simply reflect the 
workings of DM, which are, somehow, summarized 
by a very simple unifying formula? Can the ubiq- 
uitous appearance of the constant a in seemingly 
independent galactic phenomena emerge, somehow, 
from the DM paradigm? 

The Newtonian-dynamics-with-DM paradigm 



differs greatly from MOND as regards the origin and 
nature of mass discrepancies they predict in galac- 
tic systems. In MOND, these discrepancies are not 
real; the full dynamics in a given system, including 
the imaginary discrepancies, are predicted uniquely 
from the presently observed (baryonic) mass distri- 
bution. In particular, the pattern of these discrep- 
ancies is predicted, and is observed, to follow the 
well defined relations discussed in section 2. In the 
language of DM one can say that MOND predicts 
that the distribution of normal matter (baryons) 
fully determines all aspects of the DM distribution 
in each and every galactic system. In contrast, in the 
DM paradigm the expected distributions of the two 
components depend strongly on details of the par- 
ticular history of the system, since the two are sub- 
ject to different influences. The formation process of 
a given galactic system and its ensuing unknowable 
history of mergers, cannibalism, gas accretion, ejec- 
tion of baryons by supernovae and/or ram pressure 
stripping, etc., all go into determining the present 
amount, and distribution, of baryons and of DM. A 
strong and direct evidence for such differentiation 
comes from the recent realization that the ratio of 
baryon mass to the DM mass required in galaxies 
is much smaller than the cosmic value, with which 
protogalaxies should have started their way (by 
an order of magnitude, typically). This evidence 
comes, for example, from probing large galactic 
radii with weak lensing ; e.g. by Klcinhcinrich et al. 
(2004), Mandelbaum et al. (2005), and by Parker 
et al. (2007) (see also McGaugh 2007 for evidence 
based on small radius data with CDM modeling). 
Even for galaxy clusters there is now some prelimi- 
nary evidence that the observed baryon fraction is 
a few tens of percents smaller than the cosmic value 
(Afshordi et al. 2007 and references therein). 

Another acute example of the large variety of 
baryon vs. DM properties expected in the DM 
paradigm is brought into focus by the recent ob- 
servation of large mass discrepancies in three tidal- 
debris dwarf galaxies (Bournaud et al. 2007). In 
view of their specific formation scenario, hardly any 
CDM is predicted in them in the CDM paradigm. 
This is very different from what is expected in pri- 
mordial dwarfs, and contrary to what is observed. 
These dwarfs conform nicely to the predictions of 
MOND, which are based only on their presently 
observed properties (Milgrom 2007, Gentile et al. 
2007) 

How, then, can the haphazard and small amount 
of leftover baryons in galaxies determine so many of 
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the properties of the dominant DM halo, with such 
accuracy as evinced by MOND? I deem it quite 
inconceivable that DM will ever explain MOND 
and the relations it predicts for individual sys- 
tems. Achieving this within the complex scenario 
of galaxy formation in the DM paradigm would be 
akin to predicting the details of a planetary system- 
planet masses, radii, and orbits-from knowledge of 
only the properties of the central star. 

Indeed, to my knowledge none of the MOND pre- 
dictions listed in section 2 has been shown to follow 
in the DM paradigm even as just statistical correla- 
tions, to say nothing of actual predictions for indi- 
vidual systems. Even the Tully-Fisher-like relation 
that has been bruited to follow from LCDM requires 
a strong assumption on the present day baryon frac- 
tion in galaxies and on how it varies among galaxies 
(originally the ratio was assumed to be universal and 
equal to the cosmic value); but the DM paradigm 
does not come close to predicting this ratio. Not only 
isn't the general correlation predicted, but the pro- 
cesses that caused only a small fraction of the orig- 
inal baryons to show up in present day galaxies are 
likely to have produced a large scatter in this ratio, 
and hence in the predicted baryon-mass- velocity re- 
lation, unlike what is observed. 

Note also, in the context of our very subject here, 
that recent comparisons of LCDM simulation with 
observations of galaxy clusters show a clear discrep- 
ancy: observed mass distributions in the cores are 
significantly more centrally concentrated than those 
predicted by LCDM (Hennawi et al. 2007, Broad- 
hurst & Barkana 2008). 

To recapitulate, the fact that the MOND laws are 
satisfied in nature speaks against the Newtonian- 
dynamics-plus-DM paradigm in two ways: First, be- 
cause it supports MOND as a competing paradigm. 
Second, because in itself, and without reference to 
MOND, this fact directly disagrees with the expec- 
tations from the DM paradigm: A unique connection 
between the baryons and DM that holds galaxy by 
galaxy flies in the face of the haphazard baryon-DM 
relation expected in the DM paradigm. 

4. Cluster dark matter in light of MOND 

It was realized early on that MOND does not 
fully account for the mass discrepancy in galaxy 
clusters. At the advent of MOND, the baryonic 
budget of clusters was thought to consist of galaxies 
only: the large M/L values deduced for clusters-a 



few hundred (M/L) Q -were thought to represent 
mass discrepancies of order 100. The MOND cor- 
rection I found in Milgrom (1983c) reduced the 
Mj L values substantially, but still left a significant 
discrepancy (M/L of a few to a few tens (M/L) ). 
I speculated there that the x-ray emission then 
known to emanate from clusters might bespeak the 
presence of large quantities of intra-clustcr gas that 
may account for much of the remaining discrepancy. 
This has been largely borne out by the identifica- 
tion and weighing of the hot, x-ray emitting gas, 
which has reduced the cluster mass discrepancy 
(in standard physics and in MOND) by an order 
of magnitude. This, however, turned out to not 
quite suffice. Studies based on gas dynamics and 
on lensing have helped determine the remaining 
discrepancy, e.g. by The and White (1988), Gerbal 
et al. (1992), Sanders (1994,1999,2003), Aguirre et 
al. (2001), Pointccoutcau and Silk (2005), Angus 
et al. (2007b), and by Takahashi and Chiba 2007. 
In the framework of MOND we have to attribute 
the remaining discrepancy to yet undetected mat- 
ter. It is likely that this matter is baryonic in some 
form-and so I shall assume in this paper-although 
other possibilities have been raised, e.g., massive 
neutrinos (Sanders 2003, 2007). 

The following rough picture emerges regarding 
the distribution of CBDM: The observed accelera- 
tions inside a few hundred kiloparsecs of the cen- 
ter arc of the order of or a few times larger than 
dQ. Thus, MOND implies only a small correction 
there. So most of the discrepancy observed in the 
core must be due to CBDM. The ratio, I, of accumu- 
lated CBDM to visible baryons there is around 10 
to a few tens. The results of Sanders (1999) scaled 
to Hq — 75 km s _1 Mpc _1 correspond roughly to 
I = 2 at 1 Mpc. [Pointecouteau and Silk (2005) 
claim higher values. But after correcting for two 
oversights on their part their results are consistent 
with those of Sanders: for Ho = 75 km s _1 Mpc _1 
they use the small value of a appropriate for H = 
50 km s _1 Mpc _1 , and they took the mean radius in 
Sanders sample to be rather smaller than it is.] The 
value of I decrease continuously with radius, as seen, 
for example, in the small-sample study of Angus et 
al. (2007b). Since at the maximal radii in the anal- 
yses the gas mass is seen to increase faster than the 
MOND dynamical mass, we can extrapolate to even 
higher radii and conclude that for the cluster as a 
whole I is about 1 or even smaller. McGaugh (2007) 
reaches a similar conclusion based on extrapolating 
the mass- velocity relation from galaxies to clusters. 
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The contribution of the required CBDM to the to- 
tal baryonic budget in the universe is rather small. 
Fukugita and Peebles (2004) estimate the total con- 
tribution of the hot gas in clusters to to be about 
0.002, some 5 percent of the nucleosynthesis value. 
It follows from the above that the contribution of 
the CBDM is similar. 

If the CBDM is made of compact macroscopic 
objects-as is most likely-then like the galaxies, these 
must have sloshed through the hot gas for many 
dynamical times; this is implied by the CBDM dis- 
tribution being still rather extended. So, when two 
clusters collide, the CBDM will follow the galaxies 
in going through the collision zone, and will not be 
greatly affected even in head-on collisions where the 
gas components of the two clusters coalesce. Clowe 
ct al. (2006) have recently used weak lensing to map 
the mass distribution around a pair of colliding clus- 
ters and indeed found dark mass concentrations co- 
incident with the galaxy concentrations to the sides 
of the gas agglomeration near the center of collision 
(See, however, Mahdavi et al. 2007). Such observa- 
tions do not add a puzzle in the context of MOND; 
they are very much in line with what we already 
know about MOND dynamics of single clusters. An- 
gus et al. (2007a) show this with a detailed analysis. 

Another result of weak lensing analysis claimed 
to be a direct evidence for DM in clusters is the ap- 
parent ring of DM (Jee et al. 2007) in the galaxy 
cluster CI 0024+17. It turns out, however, that this 
can be explained as a MOND effect (Milgrom and 
Sanders 2007): the ring might be a direct "image" of 
the MOND transition region analogous to the GR 
transition region as marked by the Schwarzschild 
horizon. The observations and the MOND predic- 
tions are shown in Fig. 01 Note also in this connec- 
tion that gravitational lensing in MOND is more 
difficult to interpret since the thin-lens approxima- 
tion fails completely, because of the nonlinearity of 
MOND; so, the structure along the line of sight en- 
ters strongly (law 10 above). 

One may wonder why this CBDM appears in clus- 
ters but not in spiral galaxies. The unswer may lie in 
the fact that the two types of objects differ greatly 
in appearance, baryonic content, and probably in 
their formation process. Most of the observed mass 
in clusters is in the form of hot, x-ray emitting gas, 
with only a small fraction in stars, while in spirals 
the situation is the opposite. Because the required 
total amount of CBDM is comparable with that of 
the hot gas, it is conceivable that in clusters a frac- 
tion of the baryons have somehow gone into forming 



the CBDM, while this has not happened in spirals. 
If this is connected with the appearance of much of 
the baryonic content as hot gas, it may also be im- 
portant for elliptical galaxies and galaxy groups en- 
shrouded in hot gas. 

5. Heating by CBDM 

5.1. Candidates and cooling mechanisms 

Like CDM and neutrinos, some forms of CBDM 
are inherently ineffective in heating the gas. This 
would be the case, e.g., for brown dwarfs, planetary 
size objects, or too compact, cool gas clouds. On the 
other hand, black holes of different masses could ac- 
crete and radiate, or if they are massive enough they 
could lose kinetic energy to the gas by dynamical 
friction. CBDM in the form of cool dense gas clouds 
could convert their kinetic energy to thermal energy 
of the hot gas if they are large and massive enough 
to collide at a sufficient rate (e.g., Walker 1999). I 
shall concentrate on the CBDM kinetic energy as the 
source of heat: If the objects making up the CBDM 
move with the cluster virial speed, which is a reason- 
able assumption considering their distribution, the 
reservoir of their kinetic energy in the core alone is 
tens of times larger than the thermal energy of the 
gas there. If the CBDM bodies have more radial or- 
bits the core can tap even a larger reservoir. Such a 
heating source is smoothly distributed in space and 
rather steady with time. The volume density of the 
heating rate scales like the product of two densities 
(CBDM-CBDM for cloud-cloud collisions, CBDM- 
gas for dynamical friction) so there is a chance that 
the balance of the cooling rate will extend over a 
range of radii. All these are considered to be expe- 
dient for solving the puzzle (see e.g. Peterson et al. 
2003). Also, the resulting steady state temperature 
of the cooling gas has to be a fraction of the virial 
temperature and no overheating can occur. 

I now discuss briefly one possible mechanism with 
the CBDM in the form of cool dense gas clouds. 

5.2. Heating by dense gas clouds 

Cold, dense gas clouds as DM have been discussed 
by Pfenniger et al. (1994), Walker and Wardle 
(1998), and in many succeeding papers by them and 
others (see, e.g., Kamaya and Silk 2002, and ref- 
erences therein). Notably, Walker (1999) discussed 
collisions between such clouds as a dominant source 
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Fig. 4. On the left: the deduced surface density distribution around the cluster CI 0024+17 using weak lensing (Jee et al. 
2007). The small bump is the alleged ring. The other two frames show the MOND predictions for very simplistic configurations 
produced by masses well contained within the radius of the bump taken from Milgrom and Sanders (2007b). 



of the x-ray emission from galaxy clusters. Walker, 
however, assumed that these clouds constitute the 
ubiquitous DM, and also that they are responsi- 
ble for the so called extreme scattering events in 
the Galaxy. This determines the cloud parameters 
(mass and radius) quite stringently. As explained 
in the introduction, we are here totally free from 
such constraints, and can reconsider the process in 
a new light. 

Ablation of the clouds by the hot gas is also at 
work, and contributes to the heating. The rate of 
heating by ablation, or collisions, is proportional to 
the mass flux a cloud sees in his rest frame. Since 
the mass flux in CBDM is rather larger than that in 
hot gas in the cluster core (by the ratio of their den- 
sities) , one would estimate the relative contribution 
of ablation to be small. However, it can be made im- 
portant if the cross section of the clouds to ablation 
is larger than that for collisions; e.g., if clouds have 
an extended magnetic field of such strength that it 
can drag the gas but hardly affect the more compact 
clouds that fly by. In what follows I'll consider only 
collisions. 

The processes occurring today regarding the 
cloud-gas system cannot really be understood in 
disjunction from the full system history. Because, 
typically, the total CBDM mass is comparable with 
that in hot gas it is possible that, in fact, much of the 
presently observed hot gas originated in cold clouds 
that collided and converted their kinetic energy to 
heat (as proposed by Walker 1999). The evolution 
of the cloud population is thus also an important 
aspect. The life time of a cloud against destruction 
by collisions depends on its size and its location, as 
well as on properties of the cloud population. So, 



the mass and size distribution of the clouds should 
be a function of location, with much mixing on the 
cluster crossing time scale. There is also redistribu- 
tion of heat and gas on hydrodynamical time scale, 
with clouds from large radii (perhaps even from the 
very outskirts of the cluster) replenishing those at 
smaller radii where destruction is more effective. 

To get a rough idea of the properties of the clouds 
that are today most efficient in heating, I assume 
that at present all clouds are identical in size and 
mass, and require that the rate of kinetic energy 
transfer due to cloud collisions balances cooling lo- 
cally. I also assume that most of the energy released 
goes into heat, as opposed to being radiated away 
(see below). We can write the ratio of the heating 
rate, Eh, to the cooling rate, E c in a uniform spher- 
ical volume, representing the cluster core, of radius 
R = 100i? 2 kpc as: 

E H /E C ~ ii(M DM /Mg as )x 

X {TDM/Tgas - ~L){Tcool / T cr oss){^DM /S c j), (1) 

where Mjjm = 10 13 Mdm43-^o and Tdm are the 
total mass and kinetic temperature of the CBDM in 
the core volume, while M gas and T gas are those of 
the gas; r coo ; and T cross are the cooling time and the 
crossing time (= R/V, V being the 1-D velocity), 
Sum = Mdm/kR 2 is the mean surface density of 
the core CBDM, E c j = M c i/wr 2 is that of a single 
cloud (M c i being its mass, and r its radius). The co- 
efficient -q is a geometrical factor taking into account 
a Gaussian velocity distribution of the clouds, angle 
of impact, and impact parameter distribution. Using 
results from Walker (1999) I find 77 s=a 1.2. We have 
direct knowledge of all the quantities appearing on 
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the right hand side, except for £ c j, which can thus be 
constrained by the requirement for heating to bal- 
ance cooling at present. Substituting typical values 
for cooling cluster cores: r c /r cr ~ 100; Mb/M 9 ~ 
50 and {T DM /T gas ) - 3 

E d ~ 10 4 r?W ~ lO'VWia^V cm" 2 . (2) 

It is difficult to estimate the mass and size of the 
clouds separately. These may depend on how the 
clouds formed, which of them have already disap- 
pear, what supports them, etc.. 

The first rough equality in eq. @ tells us that the 
total sky cover factor of the clouds in the core is 
~ 10~ 4 , which means that the cloud-cloud-collision 
life time of a single cloud is ~ 10 4 crossing times. 
This turns out to be of the order of the Hubble time, 
which may not be a mere coincidence, but a result 
of cloud evolution that culls out clouds with shorter 
life times. The small cloud sky cover factor is also in 
line with the observations of the bullet cluster as it 
means that the CBDM masses of two colliding clus- 
ters would go through each other practically intact. 

The second rough equality in eq.([5]) tells us that 
the cloud's optical depth is very large. In a colli- 
sion, much of the kinetic energy is converted into 
thermal energy. Part of this is converted back into 
hydrodynamic energy of the gas debris expanding 
from the collision region; another part escapes as ra- 
diation. The photon escape time is typically t es ~ 
rr/c, and is at least a few times longer in our case 
than the expansion time t exp ~ r/i>, where v is the 
post-collision expansion seed, which is of the order 
of the cluster virial speed. This means that only a 
fraction of the collision energy goes into radiation 
so the energy transfer to gas heating is efficient (it 
also means that the radiation comes out at lower fre- 
quencies) . The collision debris carrying much of the 
kinetic energy of the colliding clouds will be much 
more diffused, and will readily deposit their energy 
in the gas. Calculations are needed to pinpoint the 
exact details of the collision to see if indeed all this 
can be made to work as required. 

There are, in principle, at least two ways in which 
the effects of such CBDM clouds can be observed 
directly, and which have to be studied in more de- 
tail. The first is the direct observation of radiation 
flashes coming out of individual cloud collisions. The 
characteristics of the flash: its total energy, duration, 
typical photon energy, and frequency of occurrence 
depend strongly on the mass of individual clouds 
and on the details of the collision (which, as we saw, 
involves a very optically thick system) . 



The other means of direct detection is via mi- 
crolensing of quasar images behind a cluster. Since 
the cloud surface density has to be much larger than 
lgr cm~ 2 , they are well within their Einstein radius 
for lensing at cosmological distances, and should be 
effective microlenses. The sky cover factor of Ein- 
stein discs of the clouds, which is the probability for 
appreciable microlensing amplification is, of course, 
similar to the strong lensing probability by the clus- 
ter core, which is typically not much less then unity. 
Microlensing of quasar images have been identified 
(Fohlmeister 2006a, b). The expected microlensing 
duration by clouds depends on the masses of in- 
dividual clouds. Some microlensing constraints on 
CBDM masses are mentioned by Carr (2000), but, 
as explained above should be relaxed for the CBDM 
quantities implied by MOND. 
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